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Abstract 

We study the electro-magnetic form factors of the nucleon, from small to large momentum 
transfer, in the context of the Instanton Liquid Model (ILM). As a first step, we analyze the role of 
single-instanton effects, and show that they dominate the form factors at large momentum transfer. 
Then, we go beyond the single-instanton approximation and perform a calculation to all order in the 
't Hooft interaction. We find that the ILM is in good agreement with the available experimental 
data. Based on these results, we argue that instantons provide a microscopic mechanism that 
explains the delay of the onset of the asymptotic perturbative regime, in the electro-magnetic form 
factors. 
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I. INTRODUCTION 



The recent measurements of pion and nucleon form factors performed at Jefferson Lab- 
oratory (JLAB) have triggered an important discussion about the transition from the non- 
perturbative to the perturbative regime, in QCD. The pion form factor has been measured 
very accurately from 0.6 GeV 2 < Q 2 < 1.6 GeV 2 by the F n Collaboration It was found 
that the form factor deviates significantly from the asymptotic perturbative prediction, even 
at the largest value of the momentum transfer. Important information about the proton form 
factors has been obtained by means of the recoil polarization method, which allows to access 
the ratio of the electric over magnetic form factors These experiments have shown 

that the ratio [iGe(Q 2 ) /Gm{Q 2 ) decreases very rapidly, while in the asymptotic regime it 
should approach constant 

These two results have indicated that, in elastic form factors, the asymptotic perturbative 
regime is not reached until very large values of the momentum transfer. Interestingly, this 
conclusion contrasts with the results of the CLEO experiment on 7 7* — > ttq transition form 
factor, where the asymptotic regime is reached already at Q 2 > 2 GeV 2 . For completeness, 
it should be mentioned that there exists also a combination of proton form factors which 
seems to exhibit a precocious scaling toward the perturbative behavior, namely the ratio of 
Pauli over Dirac form factors F 2 (Q 2 ) / Fi(Q 2 ) , where: 

F^Q 2 ) := ^ (G e (Q 2 ) + tG m (Q 2 )) (1) 

F2(Q 2 ) ■= (Gm(Q 2 ) - G E (Q 2 )) , r := (2) 

It was recently shown by Belitsky, Ji and Yuang that, when logarithmic corrections and 
sub- leading twist light-cone wave-functions are introduced, perturbative QCD predicts Q]: 

~ lo g 2 (Q 2 A 2 QCD )/Q 2 , (3) 

in good agreement with experimental data. On the other hand, the perturbative predictions 
for the individual Pauli and the Dirac form factor are nevertheless very far from the experi- 
mental data. This fact lead the authors to argue that the precocious scaling of the ratio (j2J) 
could be the result of a delicate cancellation in the numerator and denominator. 

The delay of the onset of the perturbative regime in the elastic form factors implies that 
there are strong non-perturbative forces inside hadrons, which dominate over the perturba- 
tive gluon-exchange even at short distances (of the order of 1/Q)- Two important theoretical 



questions arise from this fact. What is the microscopic origin of the short-scale interaction, 
driving the pion and nucleon elastic form factors away from their perturbative limit? Why 
such non-perturbative forces do not show up in the 7 7* — > tt transition form factor? 

Clearly, the answers to these questions reside in the non-perturbative sector of QCD. In 
particular, it is commonly accepted that the soft physics light quarks is very much influenced 
by the interactions responsible for chiral symmetry breaking (CSB). On the other hand, 
confinement seems to play only a marginal role. The most convincing evidence in this 
direction comes from lattice studies of QCD in the semi-classical limit: by means of the so- 
called "cooling" procedure, it was observed that, in this limit, the current-current correlators 
of light hadrons change very little, although all perturbative fluctuations are removed, and 
the string tension drops outjsj. 

The characteristic scale associated to CSB is Anf^ ~ 1.2 GeV, significantly larger than 
the typical confinement scale, Aqqd- Such a separation justifies attempting understand 
the short- distance non-perturbative structure of light-hadrons, without having simultane- 
ously to account for the microscopic origin of confinement. On the other hand, from the 
observation that 4 7r f n ~ it follows that any effective description of the short-distance 
non-perturbative dynamics of light quarks should also account for topological effects. 

Instantons are topological gauge configurations which dominate the QCD Path Integral in 
the semi-classical limit. They generate the so-called 't Hooft interaction, that solves the U(T) 

q n 

problem[7[ and spontaneously breaks chiral symmetry [8|, but does not confine. Evidence 
for such an instanton-induced interaction in QCD comes from a number of phenomenological 
studies loll , as well as from lattice simulations The ILM assumes that the 

QCD vacuum is saturated by an ensemble of instantons and anti-instantons. The only 
phenomenological parameters in the model are the average instanton size p ~ 1/3 fm and 
density n ~ 1 fm~ 4 . These values were extracted more than two decades ago, from the global 
vacuum properties [12]. The non-perturbative contribution to the electro-magnetic form 
factors of the nucleon has been analyzed in a number of works by means of phenomenological 
models (for an incomplete list see (13J). 

It the present study, we use the ILM to address the question whether the 't Hooft in- 
teraction can provide the non-perturbative dynamics needed to explain the experimental 
results on pion and nucleons's form factors. The instanton contribution to these form fac- 
tors has been investigated in the context of the ILM in a number of works. In [l4| Forkel and 
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Nielsen computed the pion form factor in a sum-rule approach, which takes into account the 
direct-instanton contribution, in addition to the lowest dimensional condensate terms in the 
operator product expansion 1 . As in other sum-rule approaches, this calculation required a 
detailed knowledge of the contribution coming from the continuum of excitations. In order 



to avoid this problem, in [16| and [17| the electro-magnetic pion and proton three-point 
functions were calculated in coordinate space, by means of numerical simulations in the 
ILM. The contribution from the continuum of excitations could be excluded by considering 
sufficiently large-sized correlation functions. The results were then compared to phenomeno- 
logical estimates of the same correlation functions, obtained from the Fourier transform of 
the fits of the experimental data. Unfortunately, this method has the shortcoming that it 
does not allow a direct comparison of the theoretical predictions against the experimental 
data. 

Direct comparison between theory and form factors at intermediate and large momentum 
transfer became possible, after the Single Instanton Approximation (SIA) was developed 

~i n n 

18, m. In 20 it was shown that instantons can quantitatively explain the pion charged 
form factor and its deviation from the perturbative regime at large momentum transfer. 
Conversely, it was observed that such effects are parametrically suppressed in the 7 7* — > 7r 
transition form factor. This explains the early onset of the perturbative regime in such a 
form factor. Moreover, a calculation of the pion distribution amplitude in the ILM was 
performed in It was found that instantons can explain the behavior of the low-energy 
experimental data (Q 2 < 2 GeV 2 ) for the 77* — > 7r° transition form factor. 

The single-instanton contribution to the nucleon electric from factors were first investi- 
gated in [22J. In this work we extend the analysis to the magnetic as well as to the Pauli 
and Dirac form factors of the nucleon. Moreover, we also go beyond the single-instanton 
approximation and include many-instanton effects, by performing a calculation to all or- 
ders in the 't Hooft interaction. We shall find that experimental data at large momentum 
transfer can be reproduced surprising well in the SIA. On the other hand, form factors at 
low momenta are very sensitive to many-instanton effects. In general, we have found very 
good agreement between theory and experiment, which indicates that instantons provide 



1 This calculation has been recently repeated, including both NLO perturbative corrections and a more 
realistic estimate of the single-instanton contribution 
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the correct non-perturbative dynamics, responsible for the electro-magnetic structure of the 
nucleon and for the delay of the onset of the perturbative asymptotic regime in elastic form 
factors. 

The paper is organized as follows. In section|n]we review the connection between the form 
factors of the nucleon and some Euclidean correlation functions, which have to be evaluated 
non-perturbatively. In section IIHI we introduce the SIA and present the predictions for 
the Sachs as well as for the Dirac and Pauli form factors of the nucleon. In section IIVI 
we include many-instanton effects by means of numerical simulations in the full-instanton 
liquid. All results are summarized in section while the appendix contains a compilation 
of the analytic SIA results. 

II. FORM FACTORS AND EUCLIDEAN CORRELATION FUNCTIONS 

In order to compute the form factors of the nucleon we consider the following Euclidean 
correlation functions: 

G3f } (t,q) = J d 3 xrf 3 ye' q - {x+y)/2 (0|Tr[^)(t,y) J 4 em (0, 0) C (n) ("^x) 74 ]|0), (4) 
G3£ re) (t,q) = J rf 3 xrfVe iq(x+y)/2 (0|Tr[^")(t,y) Jr(0,0)C (n) (-^x) 72 ]|0), (5) 

where J^ m (x) is the electro-magnetic current and r]sc n \x) is an operator which excites states 
with the quantum numbers of the nucleon. In the case of the proton we choose 2 : 

rfJx) = e abc [u T a (x)C l5 d b (x)]u c (x). (6) 

In QCD, in the limit of large Euclidean time separation t, the correlation functions (j3J and 
(J5J) relate directly to the form factors of the nucleon. In particular G3 P ^ depends linearly 
on the proton (neutron) electric form factor: 

G3t\t, q) — ► 8M 2 R(t, q) G p E {n) (Q 2 ), (7) 
R(t,q) := A 2 C e~ 2 ^, (8) 

where G^ n \Q 2 ) denotes the proton (neutron) electric form factor and A sc the coupling of 
the interpolating operator to the nucleon. Similarly, if q is chosen along the x direction, 

2 The corresponding operator for the neutron is obtained through the substitution u «-> d. 
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G3^ relates to the proton (neutron) magnetic form factor: 

G3 P # } {t, q) - -2 q 2 R(t, q) Gy"' (Q 2 ) . (9) 

These expressions are derived in the Breit-Frame, where p' = — p = q/2 and Q 2 = q 2 . 
We recall that, in the kinematic regime explored by current experiments, both Sachs form 
factors are positive definite. This implies that G3^ n '(t, q) and GS p ^ n \t, q) have opposite 
sign. 

From the correlation functions (pQ) and © it is immediate to construct linear combinations 
which relate to the Dirac and Pauli form factors: 

G3^ n) : = G3| (n) - G3?/ - 8 M 2 (1 + r) R{t, q) Ff W (Q 2 ) (10) 
G3^ n) : = -- G3S n) -G3g n) ^8M 2 (l + r) J R(t,q)F 2 p(n) (g 2 ). (11) 

Notice that, due to the sign difference, the absolute contributions of the correlation functions 
G3m and G3e to the Dirac form factor F\(Q 2 ) are added-up, while the contributions to 
the Pauli form factor ^(Q 2 ) are subtracted. 

The exponential factor R(t, q) in ©-Q and ([10 |) -([11| ) can be obtained from the two-point 
function: 

G2(t,q) = J C / 3 xe iq - x (0|Trr ?sc (t,x)r /sc (0) 7 4|0). (12) 

In the large Euclidean time limit, one has: 

G2(t,q)^2A 2 c e— \ (13) 

from which it is possible to extract the coupling constant A sc and the nucleon mass M. 

Even at asymptotically large momentum transfer, the correlation functions defined in this 
section cannot be calculated in perturbation theory. This is because all the three-point and 
two-point functions are large-sized (due to the t — >• oo limit), while pQCD is supposed to 
work only for small-sized correlation functions. On the other hand, factorization theorems 
state that, at asymptotically large values of the momentum transfer, all non-perturbative 
effects are included in the light-cone wave functions and decouple from the hard perturbative 
contributions. The problem with such an approach is that it is not possible to know a priori 
at which momenta factorization theorems become quantitatively reliable. Therefore, in this 
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work we shall refrain from using them and attempt a direct non-perturbative evaluation of 
the Green functions, from moderate to large Q 2 . 

After performing Wick contractions, the fermionically connected components of the cor- 
relators (@J) and (jHJ) read: 

G3 E(M) (i,q) = Jd 3 xJ d 3 ye iq/2 - (x+y) e a6c e oW 

ijjabca'b'c' , jjabca'b'c' , jjabca'b'c' , jjabca'b'c 1 , jjabca'b'c' , j-^abca' b' c' \ ( ~\ A\ 

\ U AA{2) + U BA{2) + U CA(2) + U D A{2) + U AA(2) + U B A{2) /) l i4 J 



where 



TJ*S$ d = Tr [S* (t, y ; -t, x) (C 75 ) T S T aa , (t, y ; -t, x) (C 75) 

S be (t, y; 0, 0)74(2)^/(0, 0; -t, x)7 4{2) ] 

^4(2?' C ' = -Tr[74(2)^(t,y;-t,x)]Tr[(C75)^ e (t,y;0,0)74(2) 

S eb , (0, 0; -t, x) (C 75 ) T S^, (t, y ; -t, x)] 
t/$r = ^[Sce(t, y; 0, 0)74(2)^(0, 0; -t, x)(C 75 ) T S cfe /(t, Y5 x)(C 75 ) T 

S L' (*> y; x ) (c 75) (t, y ; -t, x)7 4(2 )] 

= -Tr[(C 75 ) S w (f,y; -t, x)(C 75 ) T S a T a ,(t, y; -t,x)] 
Tr^^, y; 0, 0)74(2)^/(0, 0; -t, x)7 4(2 )], 



and 



Da c 4)' C ' = -Tr[ 7 4(2)^/(t,y; -t,x)] (C 75 f 5 a T a ,(t, y; -i,x)(C 75 ) 

Ste(t, y; 0, 0)74(2)5^/(0, 0; -t, x)] 
D a B b tm C ' = Tr[Sae(t, y; 0, 0) 7 4(2)S ea /(0, 0; -t, x) (C75) 
Sj b '(t,y; -t,x)7 4 ( 2 )T5f, c /(t,y; -t, x)(C7 5 ) T ]. 

Similarly, the two-point function reads: 

G2(t, p) = y d 3 xe ip - x e abc e aW < N% c a ' b ' c ' + N abc a ' b ' c ' ), (15) 

where, 

N abc a'b'c' = Tr[ 5ao/(t)X;0)0 ) (C 75) ^,(t,x;0,0)(C , 7 5 ) T ]Tr[5 cc /(t,x; 0,0)74 ] 

< c a ' bV = Tr[ S aa ,(t, x; 0, 0) (C75) S£>(f , x; 0, 0) 7 4 S b r c ,(t, x; 0, 0) {C l5 ) T ]■ 
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In these expressions, £(2/4, y; X4, x) denotes the quark propagator, the trace is over spinor 
and color indices, and the brackets ( ■ ) denote the average over all gauge field configurations. 

Fermionically disconnected components of these three-point functions bring in additional 
contribution to the form factors, coming from the quark-antiquark sea. At zero-momentum 
transfer, such contributions measure the charge of the vacuum and therefore vanish. They 
also cancel-out at finite momentum-transfer, if one assumes flavor SU(3) symmetry. 

So far, all expressions are completely general, as all the QCD dynamics resides in the quan- 
tum average over the gauge configurations. In the semi-classical limit, the non-perturbative 
contribution to the correlation functions (J3J) and © arises from single-instanton and from 
many-instanton effects. Typical single-instanton contributions are represented in Fig. ^ 
where the instanton field mediates the exchange of momentum between two partons. Many- 
instanton effects are not only those in which which a parton exchanges its momentum with 
the other partons in the nucleon by scattering on two or more pseudo-particles. In addition, 
there are also collective effects, which are associated with the breaking of chiral symmetry 
and the dynamical generation of a momentum-dependent quark effective mass jSj. These 
many-body interactions are supposed to play an important role at low momenta. 

In the next section, we shall calculate the contributions arising from the interaction of 
two massless partons with a single-instanton, while many-instanton effects will be discussed 
in section ITVl 



III. SINGLE-INSTANTON CONTRIBUTIONS 



In this section, we use the SIA to evaluate the single-instanton contribution to the cor- 
relators (UJ), © and (|12p. The SIA is an effective theory of the instanton vacuum, in which 
the degrees of freedom of the closest pseudo-particle are kept explicitly into account, while 
the contribution from all other pseudo-particles in the vacuum is included into one effective 
parameter, m* ~ 85 MeV. Such a parameter, which was rigorously defined and calculated 
in Q| for different ensembles, depends only on the two phenomenological parameters of the 
ILM, i.e. the instanton size p and density n. 

The main advantage of the SIA is that quark propagator in the single-instanton back- 
ground has a simple analytical form j^J. It consists of a zero- mode part and a non-zero 
mode part, S I (x,y) = S^. m (x,y) + S^ zm (x,y). The accuracy of the SIA was analyzed in 
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(A) (B) (C) 

FIG. 1: Graphical representation of the typical contributions to the W2W nucleon electro- magnetic 
three point function. The double lined "walls" correspond to the spatial Fourier integration. The 
dashed ellipse denotes the four quark (zero-mode) instanton interaction. The nucleon is excited at 
the left, struck by the virtual photon in the middle and absorbed at the right. Two contributions 
to the connected three-point function are shown. Diagram (A) probes the diquark content of the 
nucleon, whereas in diagram (B), the photon interacts with the remaining quark. Diagram (C) is 
disconnected, where the photon probes the sea quark content of the nucleon. 



detail in 
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181 ]. It was shown that the approach is reliable only if the relevant Green 
functions receive maximal contribution from the zero- mode part of the propagator. In fact, 
the additional 74 matrix in (jSJ) and (fT2"j) has been inserted in order to meet such a 
requirement. 

In this work, we choose to further simplify the calculation by adopting the so-called 
"zero-mode approximation", in which the non-zero mode part of the propagator is replaced 
by the free one, S T (x,y) ~ Sl^x.y) + S (x,y). Such an approximation corresponds to 
accounting for the 't Hooft interaction, and neglecting other residual instanton-induced 
interactions, which are generally sub-leading. Indeed, in Q] it was shown that the zero- 
mode approximation is very accurate in the case of the nucleon three- and two-point functions 
which we are considering. 

Finally, it is convenient to use the regular gauge and work directly in a time-momentum 
representation of the Green functions. To this end, one expresses Eqs. (114}) and ()15j) in 
terms of "wall-to-wall" (W2W) propagators, denned as the spatial Fourier transforms of the 
point-to-point (P2P) quark propagators: 

S(t',p';t,p) = J d 3 x d 3 y e' p ' y ^ p - x S(y,x). (16) 

This is achieved by insertions of appropriate delta-functions at each vertex. The convenience 



of the time-momentum representation resides in the fact that the W2 W q uark propagators 
in the single-instanton back-ground have been calculated analytically |20j and are smooth, 
non-oscillatory exponential or Bessel functions. The massless free W2W quark propagator 
is given by 

S (t', p'; t, P ) = (2vr) 3 (p' - p) 6 " um, (17) 

where w 4 = —1 and u\ = —ipi/\p\, for I = 1,2,3. The zero-mode W2W quark propagator 
in the regular gauge is given by: 

p'; t, V ) = % f(t', p'; t, p) W« A \ (18) 
m* 



f(t',p';t,p) = e *(p'-p)- K (\p'\ ~ ^) 2 + P 2 ) K (\p\ y/(t - z,f + p 2 ), (19) 
^ 7 , 7l/ i±^rjr± (20) 

where = (z, z 4 ) denotes the instanton position, m* is the effective parameter discussed 
above, and = (r, =p) are color matrices. 

The calculation of the correlation functions (@J), (0) and (|T2*|) is performed by substituting 
(I17j) and (|T8|l in the traces arising from Wick contractions. The quantum average is carried- 
out by integrating over the instanton color orientation, position, and size. The integral 
over the color orientation is trivial, while that over instanton position generates a delta- 
function which accounts for total momentum conservation. As expected, the introduction 
of an instanton-induced interaction generates an extra loop-integral, over the momentum 
exchanged through the field of the instanton. Despite the presence of loops, all diagrams 
are finite, as the instanton finite size provides a natural ultra-violet cut-off. The integral 
over the instanton size is weighted by a distribution function. In this work we assume a 
simple delta-function distribution d(p) = n5(p — p). Alternatively, one could use a fit of the 
instanton size distributions obtained from lattice simulations (for a compilation of results 
). In a previous work we have verified that these two choices essentially give the 



same result 20]. 

The SIA is reliable only if the correlation functions are dominated by the contribution of 
the closest instanton. This condition is clearly not satisfied when the distance covered by 
the quarks becomes much larger than the typical distance between two neighbor instantons. 
Previous studies [17I Q] have shown that P2P Green functions obtained analytically in the 
SIA quantitatively agree with those obtained numerically in the full instanton back-ground, 
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if the distance between the quark source and the quark sink is smaller than ~ 1 fm for 
two-point functions and than ~ 1.8 fm, for three-point functions. 

On the other hand, we do not expect the SIA calculation of the W2W correlators to be 
reliable for all values of the momentum p, even for small Euclidean times t. In fact, if the 
momentum is small the spatial Fourier transform (|l(j|) receives non-negligible contributions 
from P2P propagators connecting very distant points on the walls. However, for |p| larger 
than a GeV or so, only points at the distance smaller than roughly one inverse GeV from 
the time axis will contribute to the Fourier transform, and the SIA is applicable. 

The physical reason why at large Q 2 single-instanton effects dominate over many- 
instanton contributions is the following. In Minkowsky space, instantons correspond to 
quantum fluctuations related to tunneling between degenerate classical vacua of QCD. At 
large momentum transfer, one can imagine computing the the form factor in an infinite- 
momentum frame, where the nucleon approaches the speed of light. Following the same 
argument of Feynmann's parton model, one concludes that in this frame the dynamics of 
the nucleon is frozen. As a result of such a time dilation, during the scattering process 
quarks experience the consequences of - at most - a single tunneling event, i.e. of a single 
instanton. 

In summary, the feasibility of SIA calculations relies on the existence of a range of time 
and momentum, where the closest instanton contribution is dominant and the ground-state 
is isoMed. Previous 8tudi es UQ have shown that, for the e,etro- magn eti C three- 
point functions (jlj) and this is achieved if one chooses t to be 0.8 fm < t < 1 fm and 
restricts to the kinematic regime |p| > 1 — 2 GeV. 

In order to compare SIA predictions against experiment we shall first compute the nucleon 
coupling constant and mass from the two-point function. Then, we shall use these values to 
extract the form factors from the three-point functions. All analytic results are collected in 
the appendix. 

A. Nucleon mass and coupling constant in the SIA 

In order to extract the nucleon mass and coupling constant in the SIA, we have evaluated 
the two-point function ()12)) for t = 0.9 fm and |p| > 1 GeV. 

In Fig. Elwe show the SIA prediction for G2(t, q), and compare it with the single-pole 
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• Single-pole fit 
— SIA 
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o 1 1 1 1 1 1 1 

0.5 1 1.5 2 

|q| [GeV] 

FIG. 2: (Color online) The nucleon two-point function, with t=0.9 fm (in units of 10 6 GeV 9 ). 
The dashed line denotes the SIA prediction, the points represent the single-pole fit (|13|) . with 
A sc = 0.030 GeV 3 and M = 1.10 GeV 

fit from Eq. 0, with M — 1.10 GeV and A sc = 0.030 GeV 3 . The agreement between SIA 
and the single-pole fit is very good, indicating that, for these values of time and momenta, 
the nucleon state has been isolated. 



B. Proton form factors in the SIA 



After having extracted the nucleon mass and coupling, we are now in condition to discuss 
the single-instanton contribution to the proton form factors, which are obtained from the 
correlation functions (JZJ),©, (P"U|) and (JUJ), at 3 t = 0.9 fm. These theoretical predictions are 
affected by the errors generated by the numerical multi-dimentional loop integration and by 
the uncertainty on the best-fit values for M and A sc . The overall error is estimated to be 
smaller than 5% . 

The SIA results for the Sachs form factors of the proton are presented in Fig.EJA. and com- 



3 In (22| it was shown that, already for t > 0.7 — 1.0 fm, the relevant ratios of three- to two- point function 
are independent t. 

4 We note that the SIA prediction for Ge, reported in Fig. does not exactly coincide with the results 
reported in 22]. Such a discrepancy is due to the fact that present results are obtained with better numer- 
ical accuracy, which allowed us to determine more precisely the nucleon mass (M — 1.10 GeV ±0.01 GeV 
as opposed to the early estimate M — 1.17 GeV ± 0.05, used in [j^h 
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FIG. 3: (Color online) 

(A) SIA predictions for the proton Sachs form factors compared to experimental data 

The experimental points for electric form factor above Q 2 = 0.5 GeV 2 are obtained from the JLAB 
data for \i Ge(Q 2 )/Gm(Q 2 ), using a dipole fit for the magnetic form factor. 

(B) SIA prediction for the electric over magnetic form factors, compared to recent JLAB data 
obtained by recoil polarization method Oj. 



pared to experimental data 0, |3j, |24| . At relatively large momenta ( Q 2 > 3 GeV 2 ) , where 
the approach is supposed to work, we observe a good agreement between SIA theoretical 
calculations and experiment. 

In Fig. we show the single-instanton contribution to the the ratio of magnetic and 
electric form factors. Also in this case, we observe that theoretical calculations converge 
toward the experimental data, in the large-momentum transfer regime. However, we observe 
that at low-momentum transfer not only is the SIA curve very far from experiment, but also 
its trend is opposite. 

These results have several implications. On the one hand, we find that single-instanton 
effects provide the right amount of non-perturbative short-distance dynamics needed to 
explain the observed Sachs form factors at large momentum transfer. On the other hand, 
we see that the behavior of both the electric and the magnetic form factors at low- and 
intermediate-momentum transfer cannot be understood in terms of the interaction of the 
partons with a single instanton. In such a kinematic regime, form factors are expected to be 
very sensitive to many-instanton effects and possibly to other non-perturbative interactions. 
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FIG. 4: (Color online) 

(A) Dirac form factor of the proton evaluated in the SIA and compared to a phenomenological fit 
of the experimental data obtained as follows. In (j2j) the magnetic form factor is fitted with the 
traditional dipole formula fi" 1 = Gdi p {Q 2 ) = 1/(1 + Q 2 /0.71) 2 . The electric form factor is 
obtained from G^ t (Q 2 ) = Gdi p {Q 2 ) x (1. — 0.13 (Q 2 — 0.04)), where the second factor parametrizes 
the JLAB data for [iGm/Ge- 

(B) Q 4 times the Dirac form factor in the SIA and compared to a phenomenological fit of the 
experimental data. Perturbative QCD counting rules predict Q 4 F\(Q 2 ) ~ const. 



Let us now discuss the SIA results for the Dirac and Pauli form factors, which are re- 
ported in Fig.s (JIJ and ((HJ) and compared with the fit of the experimental results. We 
observe that single-instanton effects are sufficient to explain with impressive accuracy the 
Dirac form factor, from low- to high- Q 2 . Notice that, at the largest momentum available 
Q 2 ~ 5.6 GeV 2 , the slope of the function Q 4 F 1 (Q 2 ) is still larger than zero. On the other 
hand, we recall that in pQCD this combination should be a constant, modulo logarithmic 
corrections. Hence, we conclude that single-instanton effects provide the right amount of 
dynamics required to explain the deviation from the perturbative behavior of the Dirac form 
factor. 

The SIA prediction for the proton Pauli form factor is reported in Fig. 03 and compared 
to a fit of the experimental data. In this case, the performance of the SIA at low-momentum 
transfer is worse than in the case of the Dirac form factor. 

It is natural to ask why the same approach performs differently in the two cases. We 
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recall that the SIA is an effective theory of the ILM which can be used to account for 
instanton effects only in the limit of large momentum transfer. Therefore, the fact that the 
SIA prediction deviates from the data at small momentum tranfer does not necessarily imply 
that the instanton model is in disagreement with experiment. In order to check the ILM 
against low-energy experimental data one necessarily needs to perform a many-instanton 
calculation. 

With this in mind, let us compare the definitions of the Dirac and Pauli form factors, in 
terms of three-point correlation functions, Eq.s ()10|) and (jllj) . We observe that the F 2 {Q 2 ) 
is obtained from a difference of correlation functions of comparable magnitude (recall that 
Q3W 1 ) j g ne g a tive definite), while Fi(Q 2 ) is related to the sum of the same quantities. Notice 
also that in the combination leading to F 2 (Q 2 ) the contribution of the magnetic correlator is 
weighted by the inverse of Q 2 , (through the factor 1/r) which enhances the low-momentum 
modes, for which the SIA becomes inaccurate. From this observations it follows that the 
systematic error caused by the use of the SIA in the intermediate- and low- momentum 
regime is larger in the case of the Pauli form factor than in the case of the Dirac form factor. 

In the present calculations, all perturbative fluctuations have been neglected. It is there- 
fore important to have at least an estimate of the magnitude of these contributions. To this 
end, in Fig. we compare the complete SIA results with the predictions obtained by re- 
taining only the zero-mode part of the propagator. The difference between these two curves 
comes from free diagrams. By definition of perturbation theory, the contribution from free 
diagrams has to be larger than the perturbative corrections to them. So by comparing free 
versus zero-mode contributions, we can estimate the importance of perturbative fluctua- 
tions, relative to the non-perturbative effects we have accounted for. In all cases considered, 
we found that the instanton- induced contributions represent the dominant dynamical effect. 

In summary, we have observed that SIA is able to reproduce the Sachs form factor in 
the regime where it is applicable, i.e. at large momentum transfer. On the other hand, 
the approach misses important dynamics in the low- and intermediate-momentum transfer, 
where many-instanton effects have to be included. Interestingly, we have observed that such 
many-body contributions are not important in the Dirac form factor, which is extremely 
well reproduced in the SIA, from rather small to large Q 2 . 
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FIG. 5: (Color online) 

(A) Pauli form factor of the proton evaluated in the SIA and compared to a phenomenological fit 
of the experimental data obtained as follows. In © the magnetic form factor is fitted with the 
traditional dipole formula /i -1 Gj£ = Gdi p (Q 2 ) = 1/(1 + Q 2 /0.71) 2 . The electric form factor is 
obtained from G^ t (Q 2 ) = Gdi p {Q 2 ) x (1. — 0.13 (Q 2 — 0.04)), where the second factor parametrizes 
the JLAB data for [iGm/Ge- 

(B) Q G times the Pauli form factor of the proton in the SIA and compared to a phenomenological fit 
of the experimental data. Perturbative QCD at lowest-twist predicts Q 6 i ? 2(Q 2 ) ~ const, modulo 
logarithmic corrections. 



C. Neutron Form Factors in the SIA 



The result of the SIA calculations of the neutron electro-magnetic form factors are pre- 
sented in Fig. Eland compared with the experimental data. 

As in the case of the proton, we observe that single-instanton effects can explain the 
data on magnetic form factor in the large momentum transfer regime. On the other hand, 
the electric form factor is known only at small momentum transfer, where the SIA is not 
reliable. In this case, the SIA undershoots the experimental data by a factor two or so. 
Clearly, in order to test the validity of the ILM with such a form factor, we need to include 
many-instanton effects. 

The SIA predictions for the neutron Pauli and Dirac form factors, which are also known 
only at small momentum transfer, are presented for completeness in Fig. [7| and compared 
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FIG. 6: (Color online) 

(A) Electric form factor of the neutron evaluated in the SIA and compared to experimental 
data [23]. 

(B) Mag netic form factor of the neutron evaluated in the SIA and compared to experimental 



data 



25]. 



against experimental data. In these cases, we observe that the agreement between SIA and 
these low-energy data is indeed quite poor. 

In general, we have found that single-instanton effects alone are not sufficient to explain 
the available low-energy information on the form factors of the neutron. 



IV. MANY-INSTANTON CONTRIBUTIONS 



In the previous section we have analyzed the single-instanton contribution to the form 
factors of the nucleon. In general, we observed a good agreement with experimental data, 
in the large momentum transfer regime. On the other hand, we have verified that at low- 
momentum transfer the single-instanton effects are sub-leading, as expected. Thus, in order 
to address the question whether also the low-energy data can be explained by the 't Hooft 
interaction, we need to account for many-instanton degrees of freedom explicitly. To do so, 
we face the problem of computing the relevant correlation functions in the full instanton 
liquid vacuum, i.e. to all orders in the 't Hooft interaction. 

Such ILM calculations can be performed by exploiting the analogy between the Euclidean 



17 




Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

(A) (B) 

FIG. 7: (Color online) 

(A) Dirac form factor of the neutron evaluated in the SIA and compared to experimental data. 
The experimental curve has been obtained by assuming that the magnetic form factors follows a 
dipole formula and taking the electric form factor form experiment 

(B) Pauli form factor of the neutron evaluated in the SIA and compared to experimental data. 
The experimental curve has been obtained by assuming that the magnetic form factors follows a 
dipole formula and taking the electric form factor form experiment. 

r-i 

generating functional and the partition function of a statistical ensemble in close analogy 
with what is usually done in lattice simulations. After the integral over the fermionic degrees 
of freedom is carried out explicitly, one computes expectation values of the resulting Wick 
contractions (JUjl by performing a Montecarlo average over the configurations of an ensemble 
of instantons and anti-instantons. In the Random Instanton Liquid (RILM), the density and 
size of the pseudo-particles is kept fixed, while their position in a periodic box and their 
color orientation is generated according to a random distribution. 

In this framework, P2P correlators can be evaluated accurately in a few hours on a regu- 
lar work-station. Unfortunately, the W2W correlators which are needed in order to extract 
the form factors are much harder to compute numerically. Indeed, many simplifications 
which make the SIA approach particularly convenient do not occur in a multi-instanton 
back-ground. For example, since at the one-instanton level the W2W quark propagator 
in the instanton back-ground is known in a closed form, one can carry out calculations 
analytically, working directly in a time-momentum representation. On the other hand, in 
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a multi-instanton back-ground the quark propagator is obtained by inverting numerically 
the Dirac operator, and this operation is done in coordinate representation. Hence, one is 
left to computing numerically the six-dimensional integration in Eq.s and JHJ). Further- 
more, such an integration is complicated by the nasty oscillatory behavior of the integrand, 
introduced by the phases of the Fourier transform. 

As a result of these facts, while P2P correlators can be evaluated on an ordinary single- 
processor computer, W2W correlators typically call for a multi-processor computation. But 
even on a very powerful parallel machine, an accurate evaluation of the form factors at large 
momentum transfer is still very hard to achieve, because in such a kinematic regime the 
integrand is oscillating very fast. In this section, we propose a strategy to overcome these 
problems. 

We begin by analyzing the many-instanton contribution to the electric form factors, for 
which an important simplification occurs, as we shall see below. As a first step, we rewrite 
(J3J) as: 

G3g n) (t,q,P) = J rf 3 xrf 3 ye- 4qy+jP - x <0| Tr [tg">(2t, 0) J 4 em (*,y) C (n) (0,x) 74 ]|0). 

(21) 

Note that charge conservation implies the identity: 

G3g n) (t, 0, P) = G2(2t, P) (22) 

which can be useful to test the accuracy of the numerical integration. We can now eliminate 
one of the complex phases by setting P = 0, which corresponds to going to the nucleon's rest 
frame. At large Euclidean times, the resulting Green's function has the following spectral 
representation: 

G3g n) (t,q;0) - A s 2 c ^±^ e — ^ e - Mt G p E (n) (Q 2 ). (23) 

Now we observe that two of the three integrals in d 3 y can be performed analytically, ex- 
ploiting the fact that the above Green's function is invariant under spatial rotations 5 . We 

5 Notice that only the electric three-point function and the two-point function display such a symmetry 
property. Hence the method presented in this section cannot be applied to compute the magnetic form 
factor. 
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FIG. 8: (Color online) 

(A) : The charge distribution Green's function T3(t, |y|,0) for proton (circles) and nucleon 
(squares), evaluated in the RILM for t = 0.9 fm. 

(B) : Two-point function of the nucleon in the RILM (points), compared to a single-pole fit (dashed 
line), G2 fit (t, |q|) = 2 A 2 e -ty/i 2 + M2 with t = 0.9 fm. 



obtain: 



G3f ) (t,q,0) = i^ I %||y|sin(|q||y|)r3 p j%|y|,0) 

I ^1 1 ** 

where we have introduced the "charge distribution Green's function" : 



r3^(t,y,P) 



|^G3*%k,P) e ^. 



(24) 



(25) 



T3 p ( n \t, y, 0) represents the probability amplitude for one of the three quarks which were 
created at an initial time in a state with with quantum numbers of the proton (neutron) 
and vanishing total momentum, to absorb a photon at a distance y from the origin of the 
center of mass frame, at a later time t. When the Euclidean time becomes large such a 
Green function encodes the information about the charge distribution of the nucleon. 

Calculating numerically T3 p ^ n \t, |y|,0) for several values of |y| is not computationally 
very challenging, because it requires only a three-dimensional integration over the spatial 
position of the source and involves no oscillating phase. This problem can be handled with 
traditional adaptive Montecarlo methods, and takes a few days of computation on a regular 
single-processor machine. Then, for the final integration in d\y\, we can make use of the 
one-dimensional integration routines which are optimized for fast-oscillating functions. 
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We have evaluated the function T3(t, y, 0) in the RILM by averaging over configurations 
of 252 pseudo-particles of size p = 0.33 fm, in a periodic box 6 of volume (3.6 3 x 5.4) fm 4 . 
Like in lattice simulations, we have used a rather large current quark mass (70 MeV), in 
order to avoid finite-volume artifacts. The results for TS(t, |y|, 0) for different values of |y| 
are plotted in Fig (jSJ) A. The final one-dimensional integration in (J24|) has be handled with 
a Gauss quadrature routine, combined with a polynomial interpolation of the integrand. 

In order to extract the form factor, we have adopted the ratio of three- and two-point 



functions similar to the one suggested in 29]: 



Ge{Q ] & M + u q G2(2t,q) G2M) " (26) 
A. Nucleon Mass in the ILM 

As in the previous SIA calculation, before extracting the form factor we need to verify 
that, at the Euclidean time we work at ( t — 0.9 fm ), the contribution of the nucleon pole to 
the two-point function has been isolated. To this end, in Fig. IHlB we compare our numerical 
results in the RILM with a single-particle fit from (|13|) . The mass extracted from the fit is 
M = 1.15 GeV, in good agreement with previous estimates in the RILM 3p| . 



B. Proton Form Factors in the ILM 



The result of our calculation of the proton electric form factor in the RILM is presented 
in Fig. El where it is compared with experimental data and with the SIA curve. We observe 
a very good agreement between theory and experiment. In particular, the inclusion of many- 
instanton effects allows to explain the experimental data in the low-momentum regime, while 
at large momentum transfer the RILM gives results completely consistent with the simple 
single-instanton calculation. Quite remarkably, we find that the RILM prediction follows 
a dipole-fit at low-momenta, but falls-off faster at large momentum transfer, in agreement 
with what is observed in the recoil polarization measurements. Notice that this property 
of the form factor could not be understood at the level of the interaction of partons with a 
single-instanton (Fig. H2B). 

6 As usual, in a finite box all momenta are quantized according to pi — j- n, with n = 0, ±1, ±2, .... 
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FIG. 9: (Color online) 

(A): Electric form factor of the proton in the ILM and from experiment. Triangles are low-energy 
SLAC data, which follow a dipole fit. Circles are experimental data obtained from the recent 
JLAB result for Ge/Gm, assuming a dipole fit for the magnetic form factor. Squares are result of 
many-instanton simulations in the RILM, and the dashed line is the SIA curve. 

From the low-momentum points we can extract the proton charge radius, which falls 
slightly short of the experimental value: ( Re(rilm)) = (0-76 fm) 2 (to be compared with 
(R E r e xp \) = (0.81 fm) 2 ). The fact that we obtain a slightly small charge radius is not 
surprising. Indeed, on the one hand we recall that in the present calculation we have used 
quarks of mass of about 70 MeV, corresponding to a rather heavy nucleon (M = 1.15 GeV). 
On the other hand, we have neglected fermionically disconnected graphs, which encode 
some of the sea contribution (notice however that some "pion cloud" contribution is present 
through the Z-graphs). 

In the previous section, we have shown that the proton Dirac form factor is completely sat- 
urated by the one-instanton contribution, already at relatively low momenta (Q > 0.5 GeV 2 ). 
We can use this result, to combine the RILM result for Ge{Q 2 ) and the SIA result for 
Fi(Q 2 ) and obtain the Magnetic and the Pauli form factor of the proton in the ILM, for 
Q > 0.5 GeV 2 . The results for such form factors are reported in Fig. Also in these 
two cases we see that the inclusion of many-instanton effects is sufficient to explain the 
low-energy data. 
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FIG. 10: (Color online) 

(A) : Magnetic form factor of the proton in the ILM (squares) and from experiment (triangles). 
The ILM curve has been obtained by combining the analytical SIA prediction for the Dirac form 
factor F\ with the numerical RILM results for Ge- 

(B) : Pauli Form Factor of the proton in the ILM. The circles are obtained from a fit of the 
experimental data. Squares is the ILM prediction, obtained by combining RILM results for Ge 
with the SIA results for F\. 

C. Electric Form Factor of the Neutron in the ILM 



The results from the electric form factor of the neutron are shown in Fig. El in this 
case, the agreement with experiment is somewhat worse than the corresponding results 
for the proton electric form factor. Our theoretical prediction undershoots experimental 
data by a factor 2 or so. We believe that this discrepancy is mainly due to the absence 
of disconnected graphs. Clearly, the relative contribution of such SU(3) breaking effects 
are much more important in the case of the neutron, which has a very small electric form 
factor compared to the proton. This hypothesis is supported by the fact that in j^lj], the 
disconnected diagrams were calculated in lattice QCD and found to give a contribution of 
the order of 50% to the form factor. A systematic study of the sea contribution coming from 
disconnected graphs to several low-energy observables is currently in progress 
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FIG. 11: (Color online) 

Electric form factor of the neutron in the RILM and from experiment 2|]]. Circles are experimental 
data, squares are RILM point, while the dashed line is the SIA curve. 

V. CONCLUSIONS 

The present study was motivated by the observation that recent JLAB data show that 
electro-magnetic form factors are very sensitive to some short- distance non-perturbative 
dynamics. Instantons are known to play the leading role in the spontaneous breaking of 
chiral symmetry and in the saturation of the chiral anomaly, i.e. in two very important 
non-perturbative phenomena which occur at the GeV scale. In a previous analysis we had 
showed that instantons saturate the the pion charged form factor and, at the same time, 
explain why the perturbative regime is reached much earlier in the 77* —>■ ttq transition 
from factor. In this work, we have asked whether they can also explain some existing puzzles 
concerning the nucleon form factors. 

We have found that large momentum transfer data of Sachs as well as Pauli and Dirac 
form factors can be already reproduced by accounting for the scattering of the partons 
on a single-instanton. These calculations have been carried-out in the SIA. On the other 
hand, form factors at low-momenta cannot be calculated in the SIA because, in such a 
kinematic regime, many-instanton effects are very important. The only exception is the 
proton Dirac form factor, which is saturated by one instanton, already at relatively low 
momentum transfer (Q 2 ~ 0.5 GeV 2 ). 

We have evaluated numerically the electric form factors in the full instanton vacuum, 
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i.e. to all orders in the 't Hooft interaction, using the RILM. In the case of the proton, we 
have found that RILM predictions are consistent with SIA calculations at large momentum 
transfer and quantitatively reproduce the available body of experimental data. In particular, 
we have shown that in the ILM the electric form factor follows a dipole fit at low momenta, 
but falls-off faster at large momenta, in quantitative agreement with the recent JLAB results. 
On the other hand, the electric form factor of the neutron seems to be rather sensitive to 
fermionically disconnected graphs and SU(3) breaking effects, which have been neglected in 
the present approach. 

We have combined our SIA result for the proton Dirac form factor, with our numerical 
RILM results for the electric form factors and obtained predictions for the magnetic and 
Pauli form factors. As in the previous cases, we have found very good agreement with 
experiment for all proton form factors. In the future we are planning to use the framework 
developed in this work to investigate the role of the pion cloud in low-energy observables. 
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APPENDIX: ANALYTIC RESULTS IN THE SIA 

Every where, we choose q pointing along the 1 direction: q = (q, 0, 0). Let us define the 
following functions: 

C(t) := y/(t-ztf+(P 
e(t) := y/{t + z,f + f 

i\t) := v^T^ (27) 

A. Two-point function 

The SIA result for the two-point function defined in (jl2j) is: 

G2(t, q) = /- / dlkl dill dlml |m| 2 Ikl 2 111 2 / dc k da dc m / dz 4 (A + B), 

m* 2 7T b J J_ x 

where 

A := 2 K (\k\ r (t/2)) K (\m\ tit/2)) K Q { X ~ m ^ + (t/2)) K^x, (f/2)) e~ f < 
B := if (|k| C{t/2)) *o(M C + (t/2)) K (x;U + (*/2)) tf (x£ T(t/2)) e"**, 

and 



Xq 


■■= 


- q 2 + 2 q 1 q 


Xq 


:= >/|l|M 


- q 2 - 2 q 1 q 


Xm 


= V|l| 2 + 


m 2 + 2 m 1 \c m 




= Vlil 2 + 


m 2 — 2 m 1 c m 


xt 


= V\M 2 + 


|k| 2 + 2|k||l|c fc 


Xk ■ 


= VW + 


|k| 2 -2|k| |l| Cfc . 



B. Magnetic Three-point Function 

The SIA result for the proton (neutron) magnetic three-point function defined in (jSJ) is: 

^U M (D)(k,l, m )-±D M (U)(k,l,m) , 



G3£%q) 



n p 



rn 



*2 



d*k 



(Pi 



(2tt) 3 J (2tt) 3 J (2tt) 



m 
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where 



U M (k, 1, m) := + f/ 2 M + f/ 3 M + C/f, 

and 

£> M (k,l,m) := D? + D?. 

The functions {TfLj correspond to sets of diagrams in which the virtual photon is absorbed 
by a w-quark. They are defined as follows: 

• :=< + < + <, 



M ._ _ 256 -3^^1^^1- ^-imnm^-m ^ (| m+q | + | m |) (|k| {f)) 

|m||m + q| 

x K (\l\ K (\q/2 + m + k| £"(*)) #o(|q/2 + m + 1| 



:= +256 e ^' k+ql+2|1|) ^o(|k| ^o(|q/2 + 1 + k| r (*)) 

x K (|m| £+(*)) ^o(|q/2 + m + 1| tit)). 



< := +256 i^ii!, e -*(M+2|i|) X()( | m + q | ^ 

xX (|k| r (*)) ^o(|q/2 - 1 - k| r (*)) ^o(|q/2 + m + 1| 



11% := -1024 ,^^ 2 ^ + ^ e -*(2|q/2+m + i| + M) ^ o( | q _ k + 1 + m | r(t)) 

x K (\m + q| AT (|k| JT (|1| 



M 

2a ' ^^|q/2 + m + l||m| 



u 



_1Q24 ■ {h ~ h + ™ 2) ( ^ + h) - t (2| q /2 +m+ l| + |m +q -k + l|) ^ (|k| r(t)) 

|q/2 + m + l||q-k + m + l| Vl 1 s w; 

x X (|l + m - k| £°(f)) K Q {\\\ t(t)) X (|m| £+(*)) 
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128 

u 3a := +-, : — : n — iT (-2ZiA;i + 2l\ + 3qm 1 - 2m x k x + 4rm/i + 2m? 

Aa |m + 1 - k + q||m + q/2 + 1| v 1 1 

-qki + 2gZi + 3qh - 2l 3 k 3 + 2l\ - 2m 3 k 3 + 4m 3 l 3 + q 2 - 2l\ - 2m\ 
-Am 2 l 2 + 2m 2 k 2 + 2l 2 k 2 + 2m\ + 2 |q/2 + m + 1| |q + m + 1 - k|) 
x e - t (2| q /2 +m+ i| + | q -k + i +m |) Xo(|m| e +(i)) Xo (|l| X (|k| r (*)) *o(|m + 1 - k| 



it 



- - ~ 256 | J|| 2 q) 2 7m + 1| e ' f(2 ' m+1+q/2|+|m|) ^(|k|r(t))^o(|lk + (t)) 

X (|q + m| e(t))K (\q - k + m + 1| (0) 



|m||m + q| VM v " 

X (|q/2 + m + 1| £+(t)) K (|k| £"(*)) ^o(|q/2 + m + k| £"(*)) 



- ztM — ,.M 



m ;= _i024 v " n "* 2 - ^ ^ - ^ e- f (| m +ql+| m |)^ o( | k | r(t)) 

|m||q + m| 

X (|q/2 + m + k| r (*)) K (\q/2 + m + 1| £ + (t))^ (|l| 

The functions Df_ 2 correspond to sets of diagrams in which the photon is absorbed by a 
d-quark. They are defined as follows: 



• Df := < + d 



ib 



d M ;= -1024 , m ^ m * + l ^ i| e -*(H+2|q/2 + l + m|)^ ( | m+ [e (t)) 

|m| |m + q/2 + 1| 

K {\\\ tit)) X (|k| r (t)) K (\q - k + m + 1| C (*)) 

rf M ■= I 1Q21 — + — — - - - — -t(2|m+q/2+l|+|q+l+m-k|) 

16 ' |m + q/2 + l| + |q + l + m-k| 

K (\m + 1 - k| #o(H t(t)) K (\k\ C(t)) M\\\,e(t)) 
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jAi 



M 



jM 



< := + 256 m l e_t(|m ' +2|1|)K o(|m + q| £°(f)) *o(|k| £ + (t)) 
X (|q/2 + m + 1| e + (t)) K (\a/2 - k - 1| r (*)) 

^ := " 256 Im^^ + ilH ^'"^'^'^(Im + 1 - k + q] 

K (|m + q| *o(|k| r (0) ^o(|lU + (0) 

^ := +256 |kT^f| e ^ (lq+k|+2|1|) ^o(|l + k + q/2| r (*)) *o(H 

K (|k|e°(t))^o(|q/2 + m + lU + (t)) 

._ ^ -t(|m+l-k+q|+2|l+q/2+m|) / 1 1 1 p+ (f \\ 

d2d - 2|m + l-k + q||q/2 + m + l| e W) 
x (2m\ + 3mxg — qk\ — 2m\k\ — 2l]k\ — 2m 3 k 3 — 2l 3 k 3 + q 2 + 

3hq + 4Zimi + 2l\ + 2m^ + 4l 3 m 3 + 2/^ + 4/ 2 m 2 

-2m 2 k 2 - 2l 2 k 2 + 2l\ + 2m\ + 2|q/2 + m + 1| |q + m + 1 - k|) 

x K (\m - k + 1| £°(i)) Xo(|k| r(0) *o(H, £+(*)) 



C. Electric Three-point Function 



The proton (neutron) electric three-point function reads: 



G3g B) (f,q) 



n, 



p 4 /• d 3 k f d 3 l f d 3 



m* 2 J (2tt) 3 J (2vr) 3 J (2tt) 



m 



^^( J D)(k,l,m)-i^(f/)(k,l,m) 



where, 



U K (k, 1, m) := C/f + C/f + fXf + Uf , 



(28) 



and 



Z^(k,l,m) := Df + Di 
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As in the case of the magnetic three-point function, Uf_ A correspond to sets of diagrams in 
which the photon is absorbed by a w-quark. They are defined as follows: 

• U? :=uf a + uf b + uf c , 



< := -256 | m j| m ';y " e- X (|k| 

x K (\l\ £+(*)) K (\q/2 + m + k| C (t)) K (\q/2 + m + 1| £+(*)) 

uf b := -256 e-'d^l+^D K (\k\ X (|q/2 + 1 + k| r (*)) 

xir (|m|e + (t))^o(|q/2 + m + l|e + (t)) 



u f c : = -256 lm | +||j| 1111 e -'<H+2|H) K (\m + q| £°(f)) 
xX (|k| r (*)) ^o(|q/2 - 1 - k| r (*)) ^o(|q/2 + m + 1| 

• Ui :=ug + ug, 



ug := -1024e-'( 2 l q / 2+m+I l + l m l)K (|q-k + l + m|rW) 
x K (\m + q| *o(|k| r (0) *o(|l| 

u£ := -1024e-* (2|q/2+m+1|+|m+q - k+I|) ^ (|k|rW) 
x K (\\ + m - k| *o(|l| #o(H £+(*)) 



E ._ _| ^ff e -t(2|q/2+m+l| + |q-k+l+m|) 

3a ' |m + l-k + q||m + q/2 + l| 
x (2m 2 + 21 2 + 4ml - 21k - 2mk + 3qm 1 + 3qh + q 2 - qh 

+2 |q/2 + m + l||q + m + 1 - k|) K (\m\ 

xk (\i\ e(t)) K (\k\ r (*)) K (\m + 1 - k| e(t)) 
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ug := -256e- t ( 2 l m+1+q / 2 l + l m l)X (|k|r(0)^o(|l|e + (0) 
K (\a + m| Z°(t))K (\q - k + m + 1| C (t)) 



ug := -256 111 - q - , - 11111 111 - Hl e-^^^Kodll 
dc |m||m + q| VMS v " 

X (|q/2 + m + 1| Xo(|k| T W) *o(|q/2 + m + k| r (*)) 



(^4 •— «4a? 



;= _ 1Q24 ™ + ^ + Him + q| e - t( | m+q | +W)g (|k| 
|m||q + m| 

X (|q/2 + m + k| r (*)) K (\ q /2 + m + 1| £ + (t))^ (|l| 

The functions Df_ 2 correspond to sets of diagrams in which the photon is absorbed by a 
(i-quark. They are defined as follows: 



D? := df a + d 



ib 



df a := -1024e-' (|m|+2|q/2+I+m|) K (|m + q| K (\l\ £+(*)) 

K (\k\C(t))K (\ q -k + m + \\C(t)) 

df b := -1024 e -*(2|m + q/2+l|+| q +l+m-k|) ^ Q (| m + ! _ k | g>(fy 

K (\ m \t(t))K (\k\C(t))K (\\\e(t)) 

• iy 2 .— a 2a -i- a 2fe -|- a 2c -|- a 2d , 

4 := +256 m ^ m| l |J| m| e -«(IH+=»HI)jf (|m + q| ?(t)) K (\k\ 

K (|q/2 + m + 1| £+(t)) K (\q/2 - k - 1| r (*)) 

d| := -256 e ^(2|m+q/2 + l| + |m|)^ o( | m + J _ k + q | £- (ty 

K (\m + q| ^o(|k| r (*)) ^o(|l| (*)) 
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dg := -256 e-^+^+^Dtfofll + k + q/2| £"(*)) ^o(H £+(t)) 

Ko(|kK°(t))^o(|q/2 + m + l|e + (t)) 



2; ' () e -t(|m+l-k+q|+2|l+q/2+m|)^ o Q 1 | 



2|m + l-k + q||q/2 + m + l| 

x (2m\ + 3mig — qk\ — 2m\k\ — 2l\k\ — 2m 3 /c 3 — 2/ 3 /c 3 + q 2 + 

3hq + 4/imi + 2l\ + 2m3 + 4/ 3 m 3 + 2l\ + 4/ 2 m 2 
-2m 2 k 2 - 2l 2 k 2 + 2l\ + 2m\ + 2|q/2 + m + 1| |q + m + 1 - k|) 
x K (\m - k + 1| ^o(|k| r (*)) *o(H £+(*)) 
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